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We study numerically the influence of non-resonant effects on the dynamics of a single 



VO ! TT-pulse quantum CONTROL-NOT (CN) gate in a macroscopic ensemble of four-spin 



molecules at room temperature. The four nuclear spins in each molecule represent 
Q^ ■ a four-qubit register. The qubits are "labeled" by the characteristic frequencies, Uk, 

{k = to 3) due to the Zeeman interaction of the nuclear spins with the magnetic 
Oh! fleld. The qubits interact with each other through an Ising interaction of strength J. 

^y I The paper examines the feasibility of implementing a single-pulse quantum CN gate in 

c^ I an ensemble of quantum molecules at room temperature. We determine a parameter 

^' region, cuf^ and J, in which a single-pulse quantum CN gate can be implemented at room 

temperature. We also show that there exist characteristic critical values of parameters, 

AuJcr = \uJk' — ^k\cr and Jcr, such that for J < Jcr and AuJk = \ujk' — ujk\ < AuJcr, 
rS I non-resonant effects are sufficient to destroy the dynamics required for quantum logic 

c3 ' operations. 



I. Introduction 

Recently, the problems related to quantum computation have attracted much atten- 
tion. (See, for example, the reviews and books on this subject, [0]-[jl6|, and references 
therein.) Instead of bits with values "0" and "1" (the states of a transistor) which are 
used in a classical computer, a quantum computer uses quantum bits - qubits. Usually, 
a qubit is implemented in a quantum two-level system such as a nuclear or electron 
spin, an electron in an atom (in an ion, or in a quantum dot), or a photon in a cavity. 
The wave function of an individual qubit can be prepared in the ground state, |0), or 
in the excited state, |1), or in any superpositional state: \E' = Co|0) + Ci|l). (The only 
restriction is the normalization condition: |coP + |cip = 1.) Superpositional states 
provide significant advantages for quantum computation. 

Consider, for example, a molecule (a "register") containing four nuclear spins 
(qubits) which we shall label by "0", "1", "2", and "3", from the right to the left. 
This quantum register can be prepared initially in the state, 

^ = ^(lOa) + II3)) ® ^(102) + II2)) ® ^(|0i) + 111)) -^(|Oo) + |lo)) = 

^(lOsOaOiOo) + IO3O2O1I0) + ... + |l3l2lilo))- (1.1) 

In decimal notation, the state (1.1) can be written as, 

* = ^(|0) + |l) + ... + |15)). (1.2) 

It follows from (1.1) and (1.2) that already at the initial step of a quantum computation, 
one can load a four-qubit register with all numbers, from "0" to "15". This is not 
possible in a classical computer. 

After loading, one can perform quantum calculations using unitary transformations 
of the wave function (1.1). These calculations could utilize, for example, recently dis- 



covered Shor |17| or Grover |18] algorithm. This step of quantum computation usually 



requires error correction codes |T3| , p!9| . The final states of qubits are read using conven- 



tional techniques pOl- (See also strategies for mesoscopic measuring devices [pi|] -p7|.) 
All these operations must be completed at times smaller than the time of decoherence. 
This requirement places significant limitations on the class of materials which can be 
used for quantum computation [Q. The first quantum computation (implementation of 
Grover's algorithm) has been demonstrated using a nuclear magnetic resonance (NMR) 
techniques at room temperature [^]. In this case, the qubits were nuclear spins in a 
molecule. 

The states which are used for quantum computation are generally non-stationary 
states because they are not the eigenstates of the Hamiltonian describing the quantum 



computer. Understanding of the dynamics of these states is important. Resolving 
the dynamical issues of quantum computation is required in order to build a working 
quantum computer. One of these issues is the importance of non-resonant interactions. 
The dynamics of quantum computation involves both resonant and non-resonant inter- 
actions. Under some conditions, non-resonant effects can be ignored. In our previous 
paper |^, we studied the dynamics of the quantum CN gate in a parameter regime 



in which the non-resonant effects can be neglected. But in quantum computation 
non-resonant effects can accumulate. Even if the interaction between qubits is small, 
non-resonant effects can inhibit the desired effects and create significant errors. The 
design of a quantum computer thus requires knowing those regions of parameter space 
in which non-resonant effects are negligible. 

In this paper, we study numerically the influence of non-resonant effects on the 
dynamics of a single-pulse quantum CN gate in an ensemble of four-spin molecules 
at room temperature. We calculate the region of parameters in which a single-pulse 
quantum CN gate can be implemented at room temperature. 

The quantum CN gate is of central importance in quantum computation because 
any quantum logic gate can be decomposed into a sequence of one-qubit rotations and 
two-qubit quantum CN gates. Recent implementation of the quantum CN gate in a 
liquid state NMR involved a special sequence of resonant electromagnetic pulses. For 
example, in the "proton-carbon" nuclear spin system, the quantum CN gate was imple- 



mented using two radio-frequency pulses p9[. The first pulse induced a 7r/2 rotation 



of the target carbon spin around the one axis of the rotating reference frame. The 
second pulse induced a similar rotation around the other axis. The delay time between 
two pulses was vr/2J, where J is the interaction constant. The action of the second 
pulse depends on the state of the control proton spin. In particular, the second pulse 
cancels the action of the first pulse if the proton spin is in the ground state. In this 
experiment, the resonant frequencies were: 500 MHz for the proton spin, and 125 MHz 
for the carbon spin, J/27r = 108 Hz. In the case of a single-pulse quantum CN gate, the 
Rabi frequency, f2, must be less than the interaction constant, J, to provide a highly 
selective excitation of spin levels. For a two-pulse quantum CN gate, one should control 
the amplitude, duration and phases of two pulses, and also the delay time between the 
pulses. A single-pulse quantum CN gate is much simpler for implementation. But it 
takes more time than a two-pulse quantum CN gate because for a single-pulse quantum 
CN gate n <C J, and the pulse duration is: r = n/Vt. We show that a single-pulse 
quantum CN gate could be implemented in liquid state NMR and in solid-state sys- 
tems in which the selective excitations can be realized |1|, |^, |T2|, ^. Below we assume 



that the molecules are in a constant magnetic field. The four spins in a single molecule 
interact with each other through Ising interactions. Their nuclear spins interact also 
with electromagnetic pulses. In Section II, we present the Hamiltonian of this system, 
introduce a quantum CN gate for our system, and discuss the initial conditions for 
the density matrix. In Section HI, we present the results of numerical calculations. 



We investigate the influence of non-resonant effects on the dynamics of a quantum CN 
gate and calculate the region of parameters in which a single-pulse CN gate can be 
implemented. In the Conclusion, we summarize our results. 

II. Equations of Motion 

Consider an ensemble of non-interacting four-spin molecules which are placed in the 
magnetic field, 

B= {B^cosut,-B±smujt,B\\). (2.1) 

In (2.1), B\\ is the z-component of the magnetic field, and B± is the amplitude of the 
circularly polarized magnetic field which rotates in the {x, y) plane with the frequency, 
uj. It is convenient to represent the four spins in a single molecule by four spin indices: 
HjJ2,ki,loj where the numbers 3, 2, 1, and 0, indicate positions of nuclear spins (/ = 
1/2) in a four-qubit "register". The indices, is, J2, ki, and /q take the values, "0" 
(orientation of the corresponding spin in the positive 2;-direction), or "1" (orientation 
of the corresponding spin in the negative z-direction). All spins interact with each other 
through the Ising interaction with the constants of interaction, Ja,i3 («,/? = 0, 1,2,3). 
The Hamiltonian of the system is, 

n = H + V, (2.2) 

H = -hY: [^J: + 2 E Jc^^m]^ (2-3) 

a=0 /3>a 

h 

V = - 

In (2.2), (2.3), oja = 7a-B||, is the frequency of the spin precession due to the Zeeman 
interaction; Vt^ = 'JaBx, is the Rabi frequency; /+'~ = /^ ± il^; I^'^'^ = {l/2)a%''^'^ 
(where 7^ is the gyromagnetic ratio, and a^'^''' is the Pauli matrix, a = 0, 1, 2, 3). The 
Hamiltonian Ti is used to solve the equation for the density matrix: ihp = [7i, p] . 

Rotating system of coordinates 

The Hamiltonian 7i is time-dependent. This time dependence can be eliminated by 
transforming to a system of coordinates which rotates with the frequency of the exter- 
nal magnetic field, u. This transformation can be made using the unitary operator, 
f4; = GW[i^t(J2a=o ^a)]- ^^ ^^is rotating system of coordinates, we have the time- 
independent Hamiltonian, 

n' = ulnu^ = H' + v, (2.4) 

a=0 f3>a 




a=0 

The equation for the new density matrix, p' = UlpUuj, is: ihp' = [H',p']. To calculate 
the density matrix, p', we use the complete set of the basis states, \i3j2k1lo), 

10000) = IO3O2O1O0) = |0), 10001) = IO3O2O1I0) = |1), ..., 11111) = II3I2I1I0) = |15). 

(2.5) 
The sixteen states, \i3j2n1mQ) , form a complete set of eigenstates of the Hamiltonian 
H' in (2.4): H'\k) = El\k), {k = 0,...,15). For example, in the rotating system of 
coordinates, the energies of the ground state, |0), and the first excited state, |1), are, 

^o = -|E[(^"-^)+E'^">/3]- (2-6) 

E[ = --{-UJo + UJi + UJ2 + UJ3 - 2UJ - Jo,l - Jo,2 - Jo,3 + Jl,2 + Jl,3 + -^2,3)- 

We assume that: a;o < ^i < ^2 < ^3, and J <^ uj^+i — uj^, (« = 0, 1, 2). 

Excluding fast dynamics 

Despite the fact that the equation for the density matrix, p', has time-independent 
coefficients, this equation is still not convenient for numerical calculations. This equa- 
tion includes both fast and slow dynamics. The fast dynamics is associated with the 
presence of high frequency terms ~ (cj^ — uj). The slow dynamics is associated with 
terms proportional to the Rabi frequencies, ~ f^Q, <C cUq,. For numerical calculation 
of the density matrix, it is convenient to exclude high-frequency oscillations for each 
qubit. To do this, we use the unitary transformation: U = exp[—iH't/h). Then, we 
have the following equation for the new density matrix, 

ap" = [V",p"], {p" = U^p'U, V" = U^V'U). (2.7) 

Note that V" is a time-dependent operator. 

Initial conditions 

The thermal equilibrium initial conditions for the density matrix (2.7) are, 

PL(0)=Pfc^(0) = ^i5 ^E./k.T ^ P:U0) = 0, (^,A; = 0-15, t^k), (2.8) 

where Ek are the eigenvalues of the Hamiltonian H in (2.3) {Ek = E'^{uj = 0)). The 
expression (2.8) can be simplified using the condition, Ek/ksT <^ 1, which is usually 
satisfied for an ensemble of nuclear spins at room temperature. We assume also that 
the difference between the single-spin transition frequencies is small compared with the 



average frequency, Yl\=Q^a/^- Besides, following the idea suggested in 0], we assume 
that one can use a sequence of electromagnetic pulses which allows one to prepare the 
initial density matrix (2.8) in the state, 

p = E/16 + pA, (2.9) 




>< 0| + - - |4 >< 4| + |5 >< 5| + |6 >< 6| + |7 >< 7| + 



^^ 32kBT 

|8 >< 8| - |9 >< 9| - |10 >< 10| - 111 >< 11| J - |12 >< 12| 

In (2.9), E is the unit matrix, pA is the deviation matrix, and Tr{pj\) = 0. The idea 
suggested in 0, was to use the first four basic states, \00kl) {k, / = 0, 1), as the "active 
states" for quantum logic operations. This density matrix in the form (2.9) is chosen 
because, the dynamics of this ensemble of four-spin molecules is formally identical (up 
to some extent) to the dynamics of a pure quantum system of two interacting spins. 
Below, we shall assume (2.9) as the initial state. 

Quantum CONTROL-NOT (CN) gate 

Assume that the frequency of the external magnetic field, u, is resonant with the fre- 
quency of the transition, |0) ^-> |1), and that this transition frequency is different from 
all other single-spin transition frequencies. In this case, by applying a single vr-pulse 
with frequency uj, we shall implement the modified quantum CN gate in an ensemble 
of four-spin molecules, 

CiV = i|0000)(0001|+i|0001)(0000|+ J2 \pqrs){pqrs\= (2.10) 

p,q,r,s=0,l 
15 

i|0)(l|+i|l)(0| + ^|n)(n|, \pqrs) ^ |0000), |0001). 

n=2 

The operation (2.10) changes the state of the right-most spin (the "target" spin: "0"<-h> 
"1", only if its neighbor (the "control" spin) is in the state "0". If we assume that the 
CN gate in an ensemble of four-spin molecules works analogously to the CN gate in 
a pure quantum two-spin system, then some conditions on the density matrix must be 
satisfied. We introduce the initial wave function for the corresponding pure quantum 
two-spin system, 

\^"m = c'^mo) + c'lmi) + 4m2) + 4{om. (2.11) 

It follows from (2.10) and (2.11), that after the action of a resonant vr-pulse, one has 
the following wave function, 

l*"(r)) = c[;(r)|0) +c;'(r)|l) + 4'(r)|2) +c'3'(r)|3) = (2.12) 



CN\^"m = tc'Mll) + tc'imO) + 4'(0)|2) + 4'(0)|3), 
where the duration of a vr-pulse is approximately, 



TT 



(2.13) 



Note, that (2.13) defines the time of the vr-pulse only for a one-qubit rotation. For the 
complicated system under consideration, (2.13) gives only an approximate value for 
the time of a 7r-pulse. 

It follows from (2.12) that the following conditions should be satisfied. 



Cn(T) 



tc'liO), 



c'lir) 



<iO), 



4{r) 



c'm, 



Csl^j 



4'(o). 



(2.14) 



Now, we express the conditions (2.14) in terms of the density matrix. We introduce 
the density matrix for a pure two-spin system. 



<M^c:{0)cl"{0), (n,A; = 0-3). 



(2.15) 



/^From (2.14) and (2.15) we derive the following relation between the density matrix 
elements of the matrices, r"{T) and r"(0). 



/^ooW 



r''(r) 



r"{r) 



r'loir) 
rUr) 



01 ^ 



r''(r) 



r''(r) 



'02'. 

i2( 



r''(r) 



<Ar) 



rUr) 



rUr) 
rUr) 



20 V'; '21V'; '22V'; '23 
\rUr) ri,{r) r^r) r^r) J 



( <i(0) <o(0) ir'iM ir'i,{Q)\ 

<M r^'o(O) <2(0) <3(0) 

-<(0) -<o(0) r^'2(0) r^3(0) 

V-<(0) -^r^'o(O) r^'2(0) r^O) J 

(2.16) 

All other matrix elements of the initial density matrix, p'/^{0), in (2.9) are assumed to 

not change significantly in the process of the quantum CN operation. The relations 

(2.16) will be verified numerically in Section III for an ensemble of four-spin molecules. 

In our numerical experiment, we calculated the dynamics of the density matrix. 



PA(t) = ^i^ f: r'Ut)\n){k\, 



32kBT 



(2.17) 



n,k=0 



with the initial condition. 



^"(0)= E r„,(0)|n><A;| + if 

n,k=0 \ 



14 >< 41 + 15 >< 51 + 16 >< 61 



(2.18) 



17x71 



>< 



19x91 -llOx 101- 111 X 111 



112 X 12|. 



III. Results of Numerical Calculations 



In this section, we study numerically the dynamics of the quantum CN gate, and 
the influence of non-resonant effects on the dynamics of the quantum CN gate. We 
calculate the slow dynamics for the density matrix, r" . Below, the upper index, ", is 
omitted. The following parameters were chosen. 



cufc = cuo + 100 X A;, J„^^ = J= 10, (A;, a,/3 = 0, 1, 2, 3) 



(3.i: 



00 



{El- Eo)/h = ujo + sj, nk = n = o.i. 



(The characteristic dimensional parameters can be obtained by multiplying the values 
in (3.1) by 2n x lOS'^) The condition cu = {Ei - Eq)/Ti = Uq + SJ (or {E[ - 
Eq)/?! = 0) corresponds to the resonant transition |0) ^^ |1). As was mentioned above, 
the frequency of this transition differs from the frequencies of all other single-spin 
transitions. The density matrix in the form. 



r(0) 



/I 0\ 





Vo 0/ 



(3.2) 



corresponds to the initial state, ^E'(O) = |00), of a pure quantum system. In (3.2), we 
presented only the part of the density matrix (we shall call it below, "reduced density 
matrix" (RDM)) with the components, Vij {i,j = — 3). All other components, here 
and below, are given in (2.18). 



Dynamics of superpositional state 

The quantum CN gate (2.10) should work for any initial state, ^(0). In this sub- 
section, we consider the dynamics of the quantum CN gate (2.10) for a superpositional 
initial state. We choose here the initial RDM: 



r(0) 



/ 0.3 0.2449 

0.2449 0.2 

0.3162 0.2582 

VO.2236 0.1826 



0.3162 0.2236 \ 

0.2582 0.1826 

0.333 0.2357 

0.2357 0.1666/ 



(3.3) 



The RDM (3.3) corresponds to the following initial wave function of a pure two-spin 
quantum system. 



^(0) = Volm + v^lOl) + -^|10) + -^|11). 



(3.4) 



Now, we apply a vr-pulse with frequency, uj = ujo + 3 J. It follows from (2.16) that, 
according to the evolution of a pure quantum system, one should expect the following 



RDM at the end of the vr-pulse, 



r(t = TT/f^o) 



/ 0.2 0.2449 i0.25819 i0.1826\ 

0.2449 0.3 i0.3162 iO.2236 

-iO.2582 -i0.3162 0.333 0.2357 

V-i0.1826 -iO.2236 0.2357 0.1666/ 



(3.5) 



The RDM (3.5) corresponds to the following wave function of a pure quantum system: 
^(t = n/Vto) = iV02\00) + iV0^\01) + (1/V3)|10) + (1/V6)|11). 

The results of numerical calculations of the dynamics of the initial RDM (3.3) for 
parameters given in (3.1), are shown in Figs 1-3. As one can see from Figs 1-3, the 
values of the matrix elements of the RDM at the end of a 7r-pulse approximately coin- 
cide with the expected solution (3.5). The deviation is about 0.5%, and is caused by 
non-resonant effects. 

Influence of non-resonant effects 

As shown above, for some regions of parameters, the influence of non-resonant effects 
on the dynamics of a quantum CN gate leads only to small effects (about 0.5%). In this 
case, one can realize (at least for a finite time) quantum logic operations in an ensemble 
of four-spin molecules described by the Hamiltonian (2.3). The reason the non-resonant 
effects are small in the numerical calculations described above, is that the parameters 
(3.1) which we used, allowed us: (a) to separate significantly the frequencies, Uk, for 
different qubits (to "label" the qubits); and (b) to choose the interaction constant, J, 
between the qubits significantly large: J ^ Q. Indeed, if J <^ fi, the single-spin transi- 
tions become degenerate. For example, the frequency of the transition |0000) <-^ |0001) 
is Uo + SJ; and the frequency of the three transitions |0010) <-^ |0011), |0100) ^^ |0101), 
! ! and |1000) ^^ |1001) is Uq + J. If the value of J becomes small, all these transitions 
become degenerate, and they all contribute into the dynamics of the system. This 
degeneracy destroys quantum logic operations. Below, in this sub-section, we shall in- 
vestigate numerically the region of parameters Aujk = \ujk'—00k\ and J, which separates 
significant from insignificant infiuence of non-resonant effects on the dynamics of the 
quantum CN gate. 

The matrix elements are complex functions which are time-dependent. They include 
both time-dependent amplitudes and phases: r„fc(t) = \rnk(t)\ exp{i(f)nk(t)) , where the 
time of a vr-pulse, r, is defined in (2.13). To study the deviation of the matrix elements 
from their desired (expected) values (2.16), due to non-resonant effects, we introduce 
the deviations of the amplitudes, Aj^, and the phases, 6ij, 

Aoo(r) = ||roo(r)|-|rn(0)||, Aoi(r) = ||roi(r)|-|rio(0)||, Ao2(r) = ||ro2(r)|-|ri2(0)||, 

Ao3(r) = ||ro3(r)|-|ri3(0)||, An(r) = ||rn(r)|-|roo(0)||, A^^ir) = ||ri2(r)|-|ro2(0)||, 
Ai3(r) = ||ri3(r)|-|ro3(0)||, A22(r) = ||r22(r)|-|r22(0)||, A23(r) = ||r23(r)|-|r23(0)||. 



A33(r) = ||r33(r)|-|r33(0)||, (3.6) 

and Ajj(r) = ||rjj(r)| — |rjj(0)||, for 15 > i > 4. For the deviation of phases we have, 
according to (2.16), 

Sooir) = |0oo(r)-0ii(O)|, (5oi(r) = |0oi(r)-0io(O)|, 5o2(r) = |0o2(r)-0i2(O)-7r/2|, 

Sosir) = |0o3(r) - 0i3(O) - n/2\, 5n(r) = |0n(r) - 0oo(O)|, 
SMr) = |0i2(r) - 002(0) - 7r/2|, 5i3(r) = |0i3(r) - 0o3(O) - 7r/2|, 

S22{r) = |022(r) - 022(0)1, (523(r) = |023(r) - 023(0)1, 6s3{t) = |033(r) - 033(0)1, 

(3.7) 
and (5jj(r) = |0jj(r)|, for 15 > i,j > 4. In this sub-section the initial conditions for the 
RDM were chosen to be, 

n,{0) = n,e'-/\ r,,(0)=r*,., t < j, r,,(0) = r,„ (^,j = 0-3), (3.8) 

where ra and Vij are taken from the RDM (3.3). To study the influence of non-resonant 
effects, we introduce the dependence of the qubit frequencies, u^, on the parameter, 
M, according to the relation, 

uJk = uJo + M -k, (A; = 0,1,2,3). (3.9) 

Note, that in (3.1), M = 100, and the frequencies of qubits were well-separated. In 
Fig. 4, we show the dependence of Ajj(r) on M at the end of a 7r-pulse for the initial 
conditions (3.8), and for the following values of the parameters, 

J^^^ = J= 10, uj = {Ei-Eo)/h = uJo + 3J, Qk = ^ = 0.1. (3.10) 

The duration of the 7r-pulse is deflned in (2.13). As one can see in Fig. 4, for M > 40, 
the values of Ajj(r) are small: Ajj(r) < 0.01. As numerical calculations show, for the 
chosen values of parameters (3.10), the critical value of M is Mcr ~ 30. For M < Mcr, 
the value of Aij(r) increases as M decreases. In this case, the influence of non-resonant 
effects becomes signiflcant, and the quantum dynamics of the system described by the 
Hamiltonian (2.3) does not correspond to the quantum CN operation deflned in (2.10). 
In Fig. 5, the dependence (5ij(r) on M is shown for the same initial conditions and 
parameters used in Fig. 4. As one can see, the critical value M^r is of the same order, 
M^r ~ 30. 

In the numerical calculations presented in Figs 6 and 7, the value of M was flxed 
(M = 30) and the value of the Ising constant, J, was varied. All other parameters and 
the initial conditions were chosen as in Figs. 4 and 5. In this case, the value, M = 30 
is close to the critical value, M^r- As one can see from Figs. 6 and 7, for J > 1 the 
quantum dynamics of the system corresponds to the dynamics of the desired quantum 
CN operation within an accuracy of about 1%. For J < J^r ^ 0.3, the influence of 

10 



non-resonant effects becomes significant and one cannot use the system under consid- 
eration for quantum computation. 

4. Conclusion 

In this paper, we studied numerically the dynamics of a single-pulse quantum CN gate 
in an ensemble of four-spin molecules at room temperature. Our calculations confirm 
the equivalence between the action of a single-pulse quantum CN gate in an ensemble 
of four-spin molecules and in a pure two-spin quantum system studied earlier in |jl^ 



Chapter 25. Our results demonstrate the feasibility implementing of a single-pulse 
quantum CN gate at room temperature. 

We also have studied the influence of non-resonant effects by varying the values 
of the separation between the resonant frequencies of spins, Aw, and the constant of 
the Ising interaction, J . We have found that non-resonant effects do not destroy a 
single-pulse quantum CN gate when the ratio of Auo to the Rabi frequency, fi, is more 
than 300: /\uj/VL > 300, and for a wide range of the values of the Ising interaction 
constant, J: 3 < J/Q < 100. Consequently, a single-pulse quantum CN gate can be 
experimentally realized if the transversal relaxation time, T2, is bigger than 2n/J. 

The results of this paper can be used for experimental implementation of a single- 
pulse quantum CN gate. In particular, a single-pulse quantum CN gate could be 
implemented in a liquid state NMR at room temperature and in solid-state systems 
in which the selective excitations of the type discussed in [|I], ^ |T2|, ^ could be reahzed. 
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Figure Captions: 

Figs. 1. Dynamics of a quantum CN gate at room temperature for superpositional 
initial conditions (3.3). The matrix elements ra for i > 4 are given in (2.18). The 
parameters are presented in (3.1). A vr-pulse was applied with frequency ujq + 3 J. 
Time evolution for the first four diagonal elements of the density matrix is shown. The 
left side shows the real part of the density matrix elements. The right side shows the 
imaginary part. 

Fig. 2. Time evolution of the non-diagonal matrix elements, roi, ro2, and ro3 for 
superpositional initial conditions. The parameters and initial conditions are the same 
as in Fig. 1. 

Fig. 3. Time-evolution of the non-diagonal matrix elements, ri2, ri3, and r23 for 
superpositional initial conditions. Parameters and initial conditions are the same as in 
Fig. 1. 

Fig. 4. Influence of non-resonant effects on dynamics of a quantum CN gate. De- 
pendences of Aqj, Ai j, A2J, Aiij, and A3^3 on M, at the end of the vr-pulse. The 
initial conditions are given in (3.8). Parameters are given in (3.9) and (3.10). 

Fig. 5. Influence of non-resonant effects on dynamics of a quantum CN gate. De- 
pendences of 6j^3, and 6iij on M. Parameters and initial conditions are the same as in 
Fig. 4. 

Fig. 6. Influence of non-resonant effects on dynamics of a quantum CN gate. De- 
pendences of Aoj, Aij, A2,j, Aiij, and A3^3 on J (and M = 30), at the end of the 
TT-pulse. The initial conditions and the parameters are the same as in Fig. 4. 

Fig. 7. Influence of non-resonant effects on dynamics of a quantum CN gate. De- 
pendences of (5j^3 on J (M = 30). The initial conditions and the parameters are the 
same as in Fig. 4. 
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